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Abstract 

The low- lying excitations of the asymmetric XXZ spin chain are derived explicitly in the antifer- 
romagnetic regime through the Bethe Ansatz. It is found that a massless and conformal invariant 
phase with central charge c = 1 is separated from a massive phase by a line on which the low-lying 
excitations surprisingly scale with the lattice length as IS.E N~2, The mass gap vanishes with 
an exponent ^ as one approaches the massless phase. The connection with the asymmetric 6-vertex 
model and some physical consequences are discussed. 
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It is believed, although not rigorously proven, that the phase diagram of a classical sta- 
tistical model in d dimensions can be mapped onto that of a quantum Hamiltonian m d — 1 
dimensions In particular, singularities of the classical free energy should correspond to 
singularities of the quantum ground state energy, while the mass gap of the latter should scale, 
near criticality, as the inverse correlation length in the d-i\i direction of the classical model 
m ~ Well known examples are the 2d Ising and symmetric six- vertex models with their 
associated Ising and XXZ spin chains. In most cases, tackling the lower dimension quantum 
version turns out to be simpler, and one therefore expects the asymmetric XXZ chain to be 
more tractable than the associated asymmetric six-vertex model (i.e. with external fields). 
The six-vertex model was originally introduced to describe ferroelectric and antiferroelectric 
phase transitions in hydrogen-bonded crystals 0. Recently, there has been renewed interest 
in the solution of its asymmetric version 0, ^, due to its relation to models of crystals which 
describe the equilibrium shape of the interface between the coexisting solid and vapor phases 

The model is defined as follows: on a 2d square lattice, place on each vertical (horizontal) 
edge an up (right) or down (left) arrow. The ice condition restricts the number of allowed 
configurations to six [0], each being assigned a Boltzmann weight R^a'i^) (see fig. 1), so that 
the Yang-Baxter Equations are satisfied and the transfer matrix 

N 

forms a commuting family, [T{u),T{u')] = for any two different values of the spectral 
parameter u 0|. Introducing the Pauli matrices ± ia^) and the vertical 

polarization operator = J2k=i'^^^ ^^e associated spin chain Hamiltonian is obtained from 
the so-called extremely anisotropic limit {u 0) of the transfer matrix 

T{u) = e^^T(«) H = -log(e^^^) - sinh7;^logT(n)U=o (2) 

du 
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which gives 

^ f cosh '7 ^ 

H = - E[^^(i + ^I^I+i) + + - ^ E ^1 (3) 

where e = ±1. Since [H,T{u)] = 0, both can be diagonahzed by the same Bethe Ansatz. 
Note that for real \& (7^ 0), (|^) is nonhermitian [§. 

It is our aim in this letter to study the low-lying excitations of this asymmetric XXZ 
spin chain. Our motivation is twofold. On the one hand, the Bethe Ansatz solution for the 
asymmetric six-vertex model was originally published in a very concise letter by Sutherland, 



Yang and Yang and the phase diagram was determined from the analytical properties of 
the free energy. More details were published in 0, 0, 0, but the critical behavior has not 
been fully explored. On the other hand, the spin chain itself is of particular interest. In the 
ferromagnetic regime e = 1 and with the tuning \1/ = 7 (stochastic line), this Hamiltonian has 
been proposed as the time-evolution operator for a Id two-species diffusion process A -|- ^ 
+ y4 [1^, |15|. Under this condition H is also of significance in the understanding of Kardar- 
Parisi-Zhang-type growth phenomena |[TB[. Moreover, through a similarity transformation, 
it can be shown that (^) is equivalent to an XXZ model with boundary condition o^^_^_l = 
ex.p{±iujN)af with lu = —i^. The case of a real phase u has been studied before [0. For 
imaginary u however, the physics is completely different. 

In this letter we restrict ourselves to the antiferromagnetic regime e = —1, V = but \1' 
arbitrary. Note that V = implies taking h = —v in the six- vertex weights (fig. 1). We 
summarize our results in what follows: we solved the Bethe Ansatz in the thermodynamic 
limit and found exact, analytic expressions for the ground state energy, the mass gap and 
the dispersion relation. Our solution shows that a massive phase extends from \E' = to a 
critical "^dl) where the mass gap vanishes as m ~ (\E'c — \E')5. On the transition line the spin 
chain is gapless with a dispersion relation at small momenta AE ~ |AP|^, which indicates 
that the low-lying states scale, on a finite but large lattice, as AE ~ N^^. For \E' > \E'c 
we find a massless and conformal invariant phase, excitations scaling as AE ~ A^~^ and a 
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central charge c = 1 (Gaussian model). All these results have been checked through extensive 
numerical analysis. 

Before studying the low-lying excitations, we present some known results about the Bethe 
Ansatz equations. In a sector where n spins are flipped with respect to the reference state 
I It ... I) [T^, the Bethe Ansatz ||I9| gives for the eigenvalues of @ and their momentum 



cosh 7-2 sinh 7 ^ 0(7; a^) - V'(A^ - 2n) 

sinh 7 
cosh 7 — cos a 

j=i sinh(| + ^) 

where the rapidities {afc} are solutions of the Bethe Ansatz equations 

sinh(| + ^)]^ ^ ^^^_^^^^^sinh(7 + i(afc-«0) 



E 



0(7; a) 



(4) 



(5) 



sinh(f - ^; 



-1 



n 



(6) 



sinh(7 - - 

Eq. is found by taking the logarithmic derivative, as in (^, of eigenvalues of the transfer 
matrix. In what follows we take = 0. After taking the logarithm of (BI) we get 



1 " 27r 



with 



-% In 



/sinh(| + 



2 ' 



e(a) 



/sinh(7 + tah 
\sinh(7 — na)/ 



(7) 



(8) 



Vsinh(|-Y)y ^"Asinh(7 ^ 

{/;} is a set of n integers (half-odd numbers) if n is odd (even). Summing ^ over / and 
taking the logarithm of (13), one gets for the momentum 



P 



1=1 



N 



(9) 



The branch cuts of B(a) in the a-plane are taken from 2*7 to ioo and from —2^7 to —too 
and we choose B(0) = 0. The distribution of the integers J/ determines already a given state 
of H . It is known that at \1/ = the antiferromagnetic ground state is defined by a closely 
packed sequence of uq = N/2 integers 7° which are symmetrically distributed w.r.t. 0, and a 



set of real rapidities distributed in [— tt, tt] [|T^. As \1/ is increased, the rapidities move away 
from the real axis to form a curve C in the complex plane with endpoints —a + ib and a + ih. 
In the thermodynamic limit the integral equations determining the ground state energy are 

0(7; a)-^ f df3(j){2r, a - f3)R{f3) = R{a) (10) 



27r Jc 

lim — = Co = cosh 7 ^- — — [ daR{a)(j){'^; a) (11) 

N^oo N 2ti Jc 

Here the function R{a) is defined by 



R{ai) = hm — (12) 

N-*co N{ai+i - ai) 

It describes the density of rapidities along the curve C and is to be determined by solving (pil|). 
The endpoints of C are functions of the field \1' and the polarization y = 1 — (eigenvalue 
of jfS^). They are implicitly determined by 

p\a + ib)-^ [ d(3Q{a + ib - (3)R{(3) = -i^ + ^(1 - y) (13) 
Ztc jc I 

For \1/ less than some critical value \l/c, a = tt and —7 ( 6 ( 0, so the analytical properties 
of (|l^) allow to deform the contour C to a straight horizontal segment and the ground state 
energy remains independent of ^ [|T^] 

_ ^ ^ exp(-7|n|) 

Co = cosh 7 — smh 7 > (14) 

„±^oo cosh7n 

with the polarization remaining at y = and h given by 

* = -,-2£tH:£!f!!i^ (15) 

n cosh 7?T, 

This agrees with the fact that the free energy of the asymmetric six-vertex model does not 
depend on \1/ in this region P, [10|, |ll|. The critical value \l/c is reached at 6 = —7. Note that 



the series still converges. 

We now extend the old results to the calculation of the low-lying excitations. A complete 
study of the solutions of (||) is missing and would be desirable. Supported by numerical 



evidence we conjecture that, as for \1/ = 0, the lowest excitations are given by holes in the 



ground state distribution In particular we choose (independently of the parity of n) two 
sets of states, in the sector 77, = ^ — l=no — 1, defined by quantum numbers {//} 

h = + I (16) 

with cr = ±1, but with one If absent [Tl|] (it is understood that in the sectors n { ^ ~ 
1, one might as well consider multi-hole excitations, whose energy would be the sum of 
single hole energies. This quasi-particle feature is typical of Bethe Ansatz solvable models). 
The corresponding rapidities are shifted from their ground state set to a new set 

where one rapidity is missing (hole). To calculate the excitation energy and the difference in 



momentum between the ground state and the excited state, we define a 'shift' function [21 



jiai) = hm 4 17 

and obtain the following equations in the thermodynamic limit 

^exc -Eo = AE = -2sinh7 / da(j)'{r, «)j(«) + 2sinh7 0(7; a^'*)) (18) 

Jc 

AP = i^+p\a^''^)- [ da(f){r,a)3{a) (19) 

j{a) + ^ / dm^r, « - = ^-^&{a- a^'^) (20) 

/TT Jc I 27r 

where a^^^ is the position of the hole on the curve C and j(a;) is to be determined by solving 
equation (^. It follows from (|) and (0) that AP G [-7r,0] for a = 1 and AP G [0, vr] for 
a = -1. 

When —7 (6(0 the contour in (|T^)-(|2D|) can still be deformed to a straight horizontal 
segment, after which these equations can be dealt with in the usual way [^. When b = —7, 



to avoid the pole of ^(7, a) at a = —i-y one should close the contour by including the real 
segment [— vr, vr] and the two vertical parts along Re(a) = ±7r. We note that this alternative 
detour could be adopted for ( \l'c, too. 
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As in the case \l/ = 0, elliptic functions and elliptic integrals naturally appear in the final 
expressions. We define the elliptic modulus k by ||TT], ^ 

2 = m (21) 

TT K{k) ^ ^ 

where K(k) and K'{k) are the complete elliptic integrals of first kind with modulus k and k' 
with k"^ + /c'^ = 1 [|lT], Energy and momentum, as functions of a^^\ can be expressed in 
terms of the elliptic dn and am functions and it can be checked that Re(A£') is non negative. 
It is more instructive though to eliminate a^^"^ from the two relations (|T^) and (^) and to 
look at the single particle dispersion relation which we find to be 

A ^/ * 2 sinh7fC(A;) , , 9, , „ tto" ..l 2smh.'~iK{k) , , , 9,\1',, 1 

AE{AP) = '-^^[{1 - e sm\AP + — + —)] - + ^^[(1 - fc^ cosh2(-))] 2 

TT 2 2 TT 2 

(22) 

The mass gap is obtained by taking the minimum of this expression, i. e. by setting AP = 
(or ±7r) in (p^). This is in accordance with our numerical solution of the Bethe Ansatz 
equations which show that the lowest lying excitations are obtained when the hole in the 
distribution of rapidities is lying at the endpoints of the curve C. We obtain therefore 

m = '""^^^^'^ [(l-fe^cosh^(^))]^ (23) 

71 2 



which gives for the critical field 



cosh(y ) = k-^ (24) 



alternative to the expression obtained by taking b = —7 in (|15D. The equality of (|2|) and 
( [T3| ) for \1' = \l'c has been checked numerically. They reproduce the same critical curve. 
Fig. 2 shows this curve in the 7-\i/ plane. From (p^), for large values of 7, we obtain a 
linear dependence \1' ~ —2 In 2 + 7 + 0{e~'^), as can be seen from fig. 2. The occurence of 
the transition is intuitively clear since the terms e^a^a^ and e^*cr^cr^ tend to destroy the 
antiferromagnetic order present at = 0. 

As \E' approaches the critical value from below we obtain 

m ~ 4sinh7ir(A;)^^.^.^^^ _ ^ ^^^^ _ ^^5) 

TT 
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Alternatively, one can keep \l/ fixed and change 7, until the gap vanishes at 7c(^)- Since ( pT] ) 
defines a regular function ^(7) such that k G (0, 1) as 7 G (0, 00), we have likewise 

m ~ Co(7 - 7c(^))^ + 0(7 - 7e(^)) (26) 

The identification m ~ would therefore give a correlation length exponent = \ because 
7 is a regular function of T 0. We observe that the transition to a massless regime is 
different from the one occuring at \E' = 0, 7 = (level 1 SU (2) Wess-Zumino-Witten model 
[]23| |) where, as 7 ^ 0"^ the gap vanishes like 

m = ^ 87rexp(-7rV27) for 7 ^ 0+ (27) 

TT 

It is also instructive to compare ( p6D with the result atVy^O. From (^ it is clear that adding 
V brings a contribution —2V to the excitations in the sector n = no — 1. In this case the gap 
would vanish linearly in V at the critical value Vc = m/2, where m is given in (p^). Therefore 
(pBf ) describes a qualitatively new phase transition. 

We now study the characteristics of the model on the critical line, i. e. at \E' = \E'c. At 
small momenta we find for the dispersion relation 

4sinh7i^(A;) „ 1 , m. ,1 , ^ 

AE '-^^{k')2{l±i)\AP\2 (28) 

TT 

where positive (negative) imaginary part should be taken for AP ) (( 0). On a finite lattice, 
the momentum is quantized in units of and ( PBD suggests that low- lying excitations should 
have gaps that scale like A^~^. This behaviour has been confirmed, at least for the real part, 
by solving (^ numerically up to 80 sites and diagonalizing (|^), also numerically, up to 16 
sites. We note that the A^~^ scaling differs from the one predicted for conformally invariant 
models where the energy gaps scale as A^~^ 0|. Unusual exponents like — | for the scahng 



of the gap with the lattice length have been found in the ferromagnetic regime as well [p!4 |. 

Above the curve (p^ (see fig. 2) we find a massless and conformal invariant phase. Our 
numerical solutions show indeed that this phase is characterized by a central charge c = 1. In 
the thermodynamic limit, the ground and excited states scale as A^^^, with a compactification 
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radius i? of a free massless boson field, which depends on \l/ and 7. It is interesting to note that 
approaching the curve \1/ = "^dl) the compactification radius takes the value |. The breaking 
of conformal invariance at this value of R can probably be related to a similar phenomenon 
occurring in the XX Z model in a magnetic field {"^ = 0,V 0) where conformal invariance is 
again broken, leading to a Prokovskii-Talapov phase transition curve pS) . This phenomenon 
takes place at i? = 1 which is dual to our value i? = |. The same phenomenon occurs for 



higher spin chains |26 



Our results leads to the interesting question, namely how in the thermodynamic limit 
and in the context of the Coulomb gas picture one can account for the fact that a massless 
and conformal invariant phase ends on a curve which is also massless but exhibits different 
physics. A hint in this direction is the following observation. One expects this behavior 
to be somehow manifest in terms of a symmetry breaking. We have indeed noticed this 
numerically by diagonalizing (^) exactly for finite lattices. Above the curve (|2^) (see fig. 
2), in the massless and conformal invariant phase we find a high degeneracy in the energy 
spectrum which appears in the form of doublets and quadruplets. On the transition curve 
this degeneracy is lifted and many (but not all) of the doublets break up into singlets, while 
some quadruplets go over into a doublet plus singlets, a picture which remains valid in the 
whole massive phase. It is worth noting that this feature is already observable at finite lattice 
sizes. 

This and other points are currently being investigated. The full details of our calculations 
will be given in a future publication. 
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Figure 1: Boltzmann weights in the notation with spectral parameter u compared to that 
ref. [3]. Note that here V = j3{h + v) in contrast to the V of ref. [3] 
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Figure 2: Phase diagram in the 7-^' plane; the sohd hne corresponds to the critical line (24) 
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